Dynamical AC study of the critical behavior in Heisenberg spin glasses. 
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ABSTRACT 

We present some numerical results for the Heisenberg spin-glass model with 
Gaussian interactions, in a three dimensional cubic lattice. We measure the AC 
susceptibility as a function of temperature and determine an apparent finite tem- 
perature transition which is compatible with the chiral-glass temperature transition 
for this model. The relaxation time diverges like a power law r ~ (T — T c )~ zv 
with T c = 0.19(4) and zv = 5.0(5). Although our data indicates that the spin-glass 
transition occurs at the same temperature as the chiral glass transition, we cannot 
exclude the possibility of a chiral-spin coupling scenario for the lowest frequencies 
investigated. 
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Dynamical AC measurements in metallic spin glasses (e.g. CuMn) show the existence of 
a cusp in the in-phase component of the AC susceptibility located at a temperature value T m 
that shifts to lower temperatures as the frequency u of the AC field decreases The 
relaxation time is then given by the inverse of the AC frequency and increases very fast as 
the temperature decreases. Several laws have been proposed to describe this dependence such 
as the Vogel-Fulcher law [2j or super- Arrhenius behavior [3]. However, most experimental 
data [Sj show that T m (uS) can be well fitted to a power law over several orders of magnitude, 
T m (w) — T c ~ uj~ . T c is the value at which the characteristic relaxation time diverges and the 
power law behavior has been interpreted as a signature of a phase transition at T c . 

Most metallic spin-glasses are characterized by low spin anisotropy where a description of 
the magnetic moments in terms of continuous Heisenberg spins seems appropriate. The vast 
majority of theoretical studies in this model (analytical and numerical) have considered the 
critical behavior by studying the equilibrium properties. From these studies it emerges that 
this model undergoes a chiral-glass phase transition at a finite temperature (HI El while it was 
for long time believed that the spin-glass transition occurs only at zero temperature (Sj. A 
chiral-glass transition occurs if the spin-reflection symmetry is broken but not the spin-rotation 
symmetry. In more recent works, it was claimed that at the chiral-glass transition, there is also 
a spin-glass transition (HI E3 EEU El El • AC studies provide a direct method to investigate 
the critical behavior and are relevant as most of the experimental evidence in favour of the 
spin-glass transition is based on such type of measurements. In this letter we report on some 
dynamical AC simulations concerning the critical behavior in the Heisenberg spin-glass model 
in three dimensions. The goal of this work is then to find an estimate of the spin-glass transition 
temperature and the value of the critical exponent zv. The model is defined as 

N 

H = -J2 JijSi-Sj-h^S? , (1) 

<i,j> i=l 

where the index i runs from 1 to N = L 3 , Si is a vector of unit modulus and the < i,j > 
corresponds to a pair of nearest neighbors spins in a finite dimensional lattice with periodic 
boundary conditions. The exchange couplings Jij are taken from a Gaussian distribution with 
zero average and unit variance. Monte Carlo simulations of (^) use random updating of the 
spins with the Metropolis algorithm. A spin is randomly chosen and its value is changed as 

Oi -> <Ji + Si = ft + 5 ■ n , (2) 

with 5 a finite number (5 = 1) and fl a vector with random components extracted for a Gaussian 
distribution of unit variance. The new spin after the change in (J2J is rescaled in such a way 
that it remains of unit length. 

Our simulations of the Heisenberg model are done in the following way: An oscillating 
magnetic field h(t) = ho cos(27ru>i) of frequency u = p, where P is the period, is applied to the 
system and the magnetization measured as a function of time 

M(t) = M cos(2irut + cp) , (3) 

with Mo the intensity of the magnetization and (j) the dephasing between the magnetization 
and the field. The origin of the dephasing is dissipation in the system which prevents the 
magnetization to follow the oscillations of the magnetic field. From the magnetization we 
obtain the in-phase and out-of-phase susceptibilities defined as 

, _ M cos(^) _ 2 Jq M(t) cos(27rut)dt 

X ~ ho " h () 

M sm((p) 2 f n P M(t) sm(2Tru}t)dt . . 

X =^^ = ho ■ (5) 
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Figure 1: x' vs - T for L = 10, P = 50000 and for varying value of the amplitude of magnetic 
field h = 0.02,0.05,0.10,0.15 and 0.20 (from bottom to top). 



The dephasing cj> measures the rate of dissipation in the system and is given by 

tan(0) = ^ • (6) 
X 

The in-phase and out-of-phase susceptibilities are computed by averaging the right-hand side 
of Q, © over several periods P = —. We always averaged at least over 10 periods of time, 
after discarding the first four periods of time. We also took averages over several realizations 
of disorder, at least 6 for the largest sizes that we simulated, L = 40. 

For each frequency u, we determine the temperature T m corresponding to the maximum 
of the in-phase susceptibility x' |14| I15j . This temperature determines the relaxation time 
r(T m ) = P. Alternatively, one could also define T m as the value of the temperature at which 
an inflection point is observed in the out-of-phase susceptibility x" ■ Next, one can determine 
the static transition temperature as well as the critical exponent zv using the scaling relation 

r(T) ~ c(T - T c )- zv . (7) 

where c is a constant. In a previous work |20j . we had employed this method to study the Ising 
spin-glass model and the Heisenberg spin-glass model. Our findings, for the Ising case, was 
in very good agreement with a previous numerical study |15l I16j as well as with experimental 
results ^1E3E1EI]- For the Heisenberg case, our conclusion was that we had data compatible 
with a zero temperature divergence t(T) ~ T~ zv with a value of zv ~ 5.8. In the present work, 
we reconsider more carefully the Heisenberg case. In AC susceptibility measurements there are 
two effects that must be carefully evaluated: finite-size effects and the amplitude of the AC field. 
Finite-size effects become particularly important as we move close to the critical temperature 
where the correlation length diverges. Since one has no direct access to the correlation length, 
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Figure 2: In-phase susceptibility as function of temperature, for h = 0.01, P = 150000 and for 
L = 10, 13, 20 and 40 (from top to bottom). 



one must be sure that the size considered is large enough. The amplitude of the magnetic 
field ho is also important as we want to keep our measurements in the linear response regime. 
Already for the Ising model [TSj, it was observed that a too large value of ho can affect the 
measurements. But in the case of the Ising spin-glass model, it is only for rather large values of 
the magnetic field ho that deviations were observed in simulations, typically for ho — 0.4. Here 
on the contrary, we will see that one needs to adjust the value of ho as a function of the frequency. 
The lower the frequency, then the lower the temperature that we want to probe, and the lowest 
the magnetic field must be. For the lowest frequency that we simulated, u = 1/150000, we need 
to reduce the amplitude of the magnetic field to ho = 0.01. In Fig. ^ we show susceptibility 
values obtained for L = 10 and P = 1/uj = 50000, that is a rather small frequency. We can see 
how the position of the maximum is strongly affected by the value of ho- It is only for values 
of ho < 0.02 that the maximum converges to T ~ 0.26. However the largest value of ho up 
to which we can locate the maximum of x' f° r a given P does not seem to depend much on 
the size considered. Most of the simulation time has been spent at determining ho for each 
frequency. This was done by repeating, for each frequency, the measurements as shown in Fig^ 
for decreasing values of ho- The values of ho that we obtained are 0.05 for P < 5000, 0.02 for 
P = 15000, 50000 and 0.01 for P = 150000. 

In practice, the fact that one needs to reduce the amplitude of the applied AC field has 
an important consequence. Indeed, fluctuations in the measured susceptibility increase very 
fast as we decrease ho, as could be expected from (|4I5[) . Thus, to reduce the errors on the 
value of x'iX"> one needs to increase the number of simulated samples. Consequently, since 
ho has to be reduced as one increases/decreases the value of the period/frequency, one must 
simulate a steadily larger number of samples as the size of the system increases. In Fig. [21 we 
show x' as a function of the temperature for P = 150000, ho = 0.01 and for different sizes 
at the lowest frequency u = 1/P = 1/150000. It emerges from Fig. |2]that finite-size effects 
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Figure 3: x' vs. T for P = 1500, 5000, 15000, 50000 and 150000 for L = 40. 



are very important and strongly influence the position of the maximum of the susceptibility. 
For L = 10, the maximum is located at ~ 0.20, while for L = 13 it has moved to ~ 0.24 and 
then it stabilizes close to ~ 0.26 for L = 20, apparently not changing anymore for larger sizes. 
Thus it seems that finite-size corrections can be very important, meaning that the correlation 
length must be rather large. Although the same type of behavior is also observed at other 
frequencies, it becomes more evident as we move to lower frequencies. This might have been 
expected, as the correlation length increases when the temperature decreases. In practice, for 
all the simulations that we have performed, we always got the same susceptibility for L = 20 
and L = 40, thus we expect that our measurements for these sizes are free of finite-size effects. 

We will now present our main results for the critical behavior. We have computed the AC 
susceptibility for three sizes, L = 10, 20, 40 and for P = 1500, 5000, 15000, 50000 and 150000. 
For each of these sizes and periods, we have adjusted the value of the amplitude of the magnetic 
field ho in order to not see any shift of the maximum of x' i as explained above. Also, we have 
checked that for each value of P, the susceptibility does not change between the sizes L = 20 
and L = 40. The data that we obtained for L = 40 is shown in Fig|3J We also see how the 
value of T m decreases as u decreases. T m has been determined as the temperature for which 
there is a maximum in the susceptibility after fitting data to a parabolic form. In Fig^J we 
show the temperature T m versus the period P. Data can be fitted to a power law with 
T c = 0.19 ± 0.04 showing the existence of a finite-temperature transition in the zero-frequency 
limit. Since this value is very close to the one obtained for the chiral glass transition [7j and 
the spin-glass transition |12j . we have repeated a fit but imposing the value for T c obtained in 
these works, i.e. T c = 0.16, in order to reduce the number of parameters of the fit. With this 
condition, we obtain a value of zv = 5 ± 0.5 in good agreement with other estimates [J]. 

In this paper we have performed AC susceptibility simulations for the Heisenberg spin-glass 
model. By carefully adjusting the value of the amplitude of the applied magnetic field ho as 
a function of the frequency considered, we have determined the relaxation time associated to 
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Figure 4: Temperature of the maximum of the in-phase susceptibility vs. the period P, for 
L = 40. The discontinuous line is a best fit with (|7|). 



each temperature. Extrapolating the relaxation time to the limit of zero frequency, we extract 
a value of critical temperature T g = 0.19 ± 0.04 compatible with either the values obtained for 
the chiral-glass transition |7j or the spin-glass transition ^2 ^] . We also computed the 
exponent zv = 5.0 ± 0.5 in good agreement with previous estimates p[j and experimental data 
for Heisenberg-like models [51 1141 ET] . Thus the simplest conclusion is that we are just observing 
a spin-glass transition at the same temperature as the chiral glass transition, in agreement with 
other recent studies |3 E3- But this is not the only possibility. Indeed, in the 

coupling/decoupling spin-chirality scenario of Kawamura 0, one expects that the chirality will 
couple at small distances with the spins. Thus, as far as one considers only small distances, 
one can observe physical phenomena expressed in term of spins while the transition is really 
on the chiral parameter. It is only for distances larger than some scale L* that spins and 
chirality decouple. This crossover length is related to the value of the period or frequency 
required to probe lengths of order L* in simulations. The value of P is expected to be around 
10 5 — 10 6 |22| which is of the same order as the largest time that we probe in our simulation, 
P = 150000. Thus we cannot, with the present simulations, exclude the possibility of a chiral- 
spin coupling scenario. Simulations for much lower frequencies are needed in order to probe 
this coupling/decoupling scenario and tell if indeed, the spin transition that we observe is due 
to the chiral glass transition or not. 
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